clear all; close all; clc;

%Внутренние параметры
R1 = 220;
R2 = 220;
R3 = 220;
Rs = R1 + R2 + R3;
L = 100e-6;

%Ось времени
%fsr = L/(2*pi*Rs);
%T = 1/(10000*fsr);
%Tmod = T*10000;
T = 1e-9;
Tmod = 10^(-4);
t = 0:T:Tmod;
lt = length(t);


%Выделение памяти
Ur = zeros(1, lt);
Ul = zeros(1, lt);
i = zeros(1, lt);
di = zeros(1, lt);

%Внешние параметры, внешние фазовые переменные
%E = 10*ones(1, lt);
E = 0*(t<Tmod/4) + 3*(t>=Tmod/4);

%/Начальное состояние
i(1) = 0;
di(1) = 0;
Ul(1) = di(1)*L;

for k = 2:lt
    Ul(k) = E(k) - i(k-1)*Rs;
    di(k) = (Ul(k))/L;
    i(k) = i(k-1)+di(k-1)*T;
    Ur(k) = i(k)*Rs;
end

figure(1);
plot(t*1e9, [E;Ur]);
xlabel('t, ns')
ylabel('E, U, Volt')
legend('U(t)', 'E(t)')
grid on

